Exotic ground states and impurities in multiband superconductors 
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^\ i We consider the effect of isotropic impurity scattering on the exotic superconducting states 

t-H . that arise from the usual BCS mechanism in substances of cubic and hexagonal symmetry where 

i-H ' the Fermi surface contains inequivalent but degenerate pockets {e.g. around several points of high 

symmetry). As examples we look at CeCo2, CeRu2, and LaBe; all of which have such Fermi surface 
topologies and the former exhibits unconventional superconducting behavior. We find that while 
these non s-wave states are suppressed by non-magnetic impurities, the suppression is much weaker 
CI than would be expected for unconventional superconductors with isotropic non-magnetic impurity 

scattering. 

Recently it was shown that in substances with cubic or hexagonal symmetry, exotic superconducting states can 
O ' arise from conventional BCS mechanisms This can occur when these metals have Fermi surface (FS) pockets 
^ | that are centered at or around some three- fold or higher degenerate symmetry points of the Brillouin zone (BZ). 
5h ' The resulting superconducting state is either a conventional s-wave state or corresponds to a multidimensional su- 
Q*', perconducting representation for which the ground state breaks time reversal symmetry. Such FS topologies exist 
' in many superconductors. Examples are CeRu 2 0, LaB 6 j|, and CeCo 2 0; the latter exhibits a low temperature 
, t | behavior specific heat that appears to vary as a power law with temperature jj^] which is consistent with unconven- 
3 ■ tional superconductivity. Here we investigate the stability of these exotic superconducting states in the presence of 
non-magnetic isotropic impurity scattering. In particular we examine in detail a model that applies both to CeCo2 
and to the FS pockets centered at the six N-points of a body-centered cubic (BCC) lattice and discuss the results of 
the analogous theory for CeRu2 and LaB6- We show that the exotic superconducting states are less susceptible to 
non-magnetic impurities than might be expected from earlier theories of isotropic impurity scattering in single band 
O | unconventional superconductors (|-@] • 

A BCS approximation for the multi-band case will be used (see e.g. ||). The form of the interaction parameters 
T—i \ describing the two electron scattering on and between the different degenerate pockets of the FS is fixed by symmetry. 
^ • Consequently, the resulting superconducting state need not be s-wave when the interaction does not depend on the 
direction of the Fermi momenta of an individual pocket. In our earlier work Q we considered three cases in detail a) 
three FS pockets centered about the X-points of a simple cubic lattice; b) three FS pockets at the M-points of the 

j. hexagonal lattice; c) four FS pockets at the L-points in the face-centered cubic (FCC) lattice. Case a) may describe 

the superconducting states in LaB 6 and in CeRu2 since these materials both have FS pockets centered at the X-point. 
0^ ' It has recently been demonstrated that the specific heat exhibits a power law temperature dependence at low 
temperatures in CeCo2 @||. For this reason we consider here FS pockets centered at the six N-points of the BCC 
lattice which is formally identical to that in FCC CeCo2 for the FS pockets which are located along the (1,1,0) and 
equivalent directions [|||] (note that in this material there also exist FS sheets centered at the T point). 
The Hamiltonian for several separate pieces of the FS can be written in the following form: 

"6 1 1 
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where a and a' are spin indices, A Q ^(q) includes the interaction for scattering two electrons from the pocket a 
into pocket (3 which is due to both Coulomb and electron-phonon terms. For simplicity we will take A Q ./3(q) to be 
independent of the direction of q. Scattering by impurities is described by the Hamiltonian 

[i(q - k)R J -]V a , /3 (k,q)aJ (r (k) £ w(q). (2) 

(7 k,q a,/3 j 

Introducing the anomalous Green's function J- a (x — x') for each FS sheet a, the corresponding Gor'kov equations can 
be found for the case of singlet pairing in exactly the same manner as in Ref [ [To| . Including the isotropic impurity 
scattering within the Born approximation the following Gor'kov equations are found 

[iu> n - e Q (k) - G a (oj)}g a (k,iuj n ) - [A* + ^l(w)]J r Q (k, iu n ) = 1 
[iw„ + e a (k) + Cj a (-u;)]J-t(k, iu n ) - [A* + J^(w)]a a (k, iw„) = 



1 



where 



G a (u) = ^E \u a ,p\ 2 Gp{Kiu n ) 



k,/3 n 



\ u a,p\ 2 = J J ~r~\Vot,pO s -Fi k^)| 2 and rii is the concentration of impurities. Introducing A Qi „ = A a + fl(w) and 
Wa,n = iw n ~ the gap and self-energy equations can be expressed as 
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where r Q / 3 = 7r7Vo|it Q ./3| 2 , and Nq is the normal density of states on a single pocket. Equation s (3|, |i| ), and @ have 
been studied in the context of multi-band superconductivity by a variety of authors (see, e.g., 1 1 - 1 3| ) . 

We consider here FS pockets centered at the six N points of the BCC lattice. This situation is formally identical to 
that describing the FS pockets in FCC CeCo2 located along the (110) and equivalent directions. The N points of the 
BCC lattice lie at ib$ and i(bj — hj) with i > j where the are the reciprocal lattice basis vectors [bi = ^(Oj 1j 1)> 
l>2 = —(1, 0, 1), and b3 = — (1, 1, 0)]. The interaction between the electrons forming a Cooper pair on Fermi pockets 
centered at these points takes the form 
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where A is the interaction on the same pocket, /i couples two pockets connected by a —(1,0,0), —(0,1,0), or a 
-^-(0,0,1) translation, and v characterizes the remaining couplings between nearest neighbor N points. The matrix 
for the impurity scattering (with elements r Qij g) has the same form and will be described by elements ro,ri, and T2 
replacing A, is, and fi respectively. 

Consider the clean limit, the linearized gap equation is 



(J) 



The six A a transform among each other under cubic symmetry transformations forming a 6D reducible representation 
of the cubic group Oft, which is split into a ID A lg , a 2D E g , and a 3D F 2g irreducible representation. These three 
representations correspond to different order parameters with three critical temperatures: 



J c o,A = exp 



T, 



cO,E 



T C 0,F — 



7T 

2-fUJp 
7T 



1 



exp 
exp 



N a (\ + n + 4is) 
1 

N (X + n-2iy) 
1 



(8) 
(9) 
(10) 



2 



where the factors in the exponentials must be negative for a non-zero transition temperature. When v > or 
2v + fi > {i.e. for repulsive inter-pocket interactions) the higher dimensional representations have the higher T c . 
The basis wave function for ID A\ g identical representation is 

I = (Ax + A 2 + A 3 + A 4 + A 5 + A 6 )/V6, (11) 
the basis wave functions for the 2D E g representation can be chosen as 

m = (Ai + eA 2 + A 3 + eA 4 + e 2 A 5 + e 2 A 6 )/%/6 

r/ 2 = (Ai + e 2 A 2 + A 3 + e 2 A 4 + eA 5 + eA 6 )/V6 (12) 
where e = exp(27ri/3), and the basis wave functions for the 3D F% g representation can be chosen as 

7} x = (A s - A 6 )/V2 

Vy = (A x - A 3 )/V2 (13) 
7] Z = (A 2 - A 4 )/V2. 

Following Refs. Jl4|,[l|], the Ginzburg Landau functional for the 2D and 3D representations are found to be 

6f e /n = T ~ Tc0 ' E (\vx\ 2 + \V2\ 2 ) + Jif 2 (1ml 4 + \m\ A + MmW) (14) 

1 cQ : E 9D7T z i c0E 



SF F /N = T T Tc °' F (rfrT) + "™ (fe| 4 + k| 4 + fe| 4 )- (15) 
J-cO,F 62tt z 1 cQF 

Eq. ( fl4|) for the E g representation implies (771,772) = (1, 0) is a stable ground state. Following the notation of Ref. fl5|] 
this corresponds to the superconducting class 0(1)2) ■ The properties of such a ground state are the same as those 
found for the analogous calculation for FS pockets centered at the three X points in a simple cubic lattice and this is 
also the situation that applies to both LaBg and CeRu2 (we refer to Ref. Q for more details). Eq. ([l5|) for the F 2g 
representation implies the ground state solution (r] Xj r]y,r] z ) = [1, exp(i<p±), exp(i^>2)] where tfti and 4>2 are arbitrary 
phase factors. The degeneracy arising from the phases 4>\ and $2 is an artifact of the BCS theory and is not lifted 
by higher order terms in free energy functional. In the notation of Ref. |Ie| the free energy of Eq. (jl^) places the 
system on boundary of two phases with superconducting classes D3 x R and D^{E). A similar situation arose for FS 
pockets centered at the four L points of a face centered cubic lattice where a degeneracy was found between solutions 
corresponding to the classes (D2) x R and D 4 (i?) of the F2 g representation 0. The presence of a FS centered at 
the r point lifts this degeneracy. As a result the solution {rj x ,rj yi rj z ) — (1, e, e 2 ) corresponding to the magnetic class 
Ds(E) is likely. This class allows ferromagnetism and it will have point nodes on the T-centered FS (the N point 
centered FS pockets will have no nodes) |l^,|l6| . 

Now consider the effect of including impurity scattering on T c ^AtT c .f,T c ,e- It is convenient to define the following 
relaxation times 

(2r A r X = To + 4Ti + T 2 

{2TE)- 1 = r + r 2 -2r 1 (16) 
(2t f )- x = r - r 2 (17) 

The lifetime ta determines the elastic mean free path while the lifetimes te and t f cannot be easily measured. 
Solution of the linearized equations indicate that the transition temperatures T C) j are given by 

b l7 = ^ + I w ) -^ ) (18) 

where ip(x) is the digamma function. From this expression it is clear that T c a is not suppressed by non-magnetic 
impurities while T c .e and T c _p are suppressed. Note that if T± = T2 = (there is only intra-pocket impurity scattering) 
then all three transition temperatures are not suppressed by impurities. It is only inter-pocket scattering that reduce 
T Cj e and T Cj f from T c o,b and T c o,_f. Consequently, the suppression of T c as the impurity concentration is increased 
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is not as rapid as might be expected from previous analysis of unconventional single band superconductors HQ. 
Furthermore an impurity induced transition is possible from the non s-wave superconducting states to the s-wave 
superconducting states when A < —Av — fx where negative A corresponds to an attractive intra-pocket interaction and 
either v > or 2v + fi > 0. The observation of an abrupt change in the T c dependence on impurity concentration 
corresponding to this phase transition (i.e. T c initially decreasing with impurity concentration in the exotic state and 
then becoming impurity independent in the s-wave state) may help to identify such exotic superconducting states. 

The results of the above analysis may be directly applied to CeCo2 and predicts that due to the presence of 
the the T-centered FS sheet, the low temperature specific heat will vary as T 3 for the exotic states considered (with 
superconducting classes 0(D 2 ) or D 3 (E)). Note that initial measurements appear more consistent with a T 2 behavior 
however these measurements only exist down to T/T c 0.13 Q which may not yet be sufficiently low to extract the 
low temperature behavior reliably. For the materials LaBg and CeRu2 the FS pockets are centered at the X points for 
which the exotic superconducting state belongs to the 2D E g representation which is described in detail in Ref. [Q. 
Impurity scattering will have qualitatively the same effect on this state as on the exotic states examined above. 

In conclusion we have considered the effect of isotropic non-magnetic impurity scattering on the exotic supercon- 
ducting states that arise in a BCC lattice when the FS forms pockets centered at the six N-points. We have found 
that increasing the impurity concentration does not suppress these non s-wave states as rapidly as would be expected 
from earlier studies of impurity scattering in unconventional superconductors. These results are not specific to FS 
pockets centered at the reciprocal lattice N points of a BCC lattice but will apply to any degenerate set of FS pockets. 

We would like to thank Z.Fisk, D. Khokhlov, J.R. Schrieffer, and the members of the NHMFL condensed matter 
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Field Laboratory through NSF cooperative agreement No. DMR-9527035 and the State of Florida. 
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